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"NITRIC  OXIDE  PRODUCTION  IN  LAMINAR  AND  TURBULENT 
DIFFUSION  FLAMES" 


Arnable  Lilian* 

ABSTRACT 

An  analysis  of  the  production  of  Nitric  Oxide  in  laminar 
and  turbulent  diffusion  flames  has  been  carried  out  taking 
advantage  of  the  large  activation  energy  of  the  NO  producing 
reactions . 

In  laminar  diffusion  flames  the  NO  production  takes  place 
in  a thin  high  temperature  region  adjacent  to  the.  thin  reaction 
zone  for  the  main  reaction.  The  mass  of  NO  produced  per  unit 
flame  surface  and  time  is  calculated  in  terms  of  the  main  charac- 
teristics of  the  Burke-Schuinann  thin  flame  solution. 

The  NO  production  rate  in  turbulent  diffusion  flames  is 
related  to  the  probability  density  function  for  an  inert  species 


Universidad  Politecnieu  de  Madrid,  Spain. 


A 


I 


INTRODUCTION 


The  production  of  Nitric  Oxide  has  been  analyzed  for  the 
laminar  diffusion  flames  associated  with  fuel  droplet  combustion 
by  Bracco  ( 1 972  ),  Altcnkirch  et  ul . (1972  ) and  Kesten  ( 1972  ). 

The  "HO"  production  in  turbulent  diffusion  flames  has  been  anal- 
yzed, among  others,  by  Fletcher  and  Hcywood  (1971),  Quang  et  a.l  . 
(1  972  ),  Kent  and  Bilger  ( 1972  ) and  (1  974  ) and  by  Gouldin  (197*0. 

In  most  of  these  studies  the  instantaneous  rate  of  "NO" 
production  has  been  determined  by  using  the  reaction  mechanism 
proposed  by  Zeldovich  (1946) 


(0)  + 

K l 

(N?)  -►  (NO)  + 

(N) 

K?. 

Cl)  + 

( 0 2 ) ->•  (NO)  + 

(0) 

If  the  concentration  of  NO  in  the  region  of  IT 0 production  is 
small  compared  with  its  local  equilibrium  value  the  effects  of 
the  reverse  reactions  can  be  neglected.  The  rate  constant  fev 
the  second  reaction  is  large  compared  with  that  of  the  first 
reaction  so  that  the  nitrogen  atom  concentration  N can  be  obtained 
using  the  steady  state  approximation 


k^Ng)  {0}  = k2(N>  { 02 } (3) 

in  terms  of  the  molar  concentrations  of  N^,  0 and  0,  and,  in 
addition,  the  production  rate  of  NO  in  moles/cm  sec.  is  given 
by 

{NO}  = 2k1  { N 2 } {0}  (3a) 

3 

where  according  to  Wray  and  Tear e (3 962  ) , is  given  in  cm  /moJe.sec 


7.10  exp  (-38100/T  ) , with  T in 


by  k1 


o 1/ 


The  rate  constant  k ^ is  very  often  small  compared  with 
those  of  the  reactions  between  the  main  reacting  species  so  1 hat 
the  changes  in  concentrations  of  Ng  > 0^  an<^  0 associated  with 
the  reactions  (1)  and  (2)  can  be  neglected  coinjKtrcd  those  associ- 
ated with  the  main  reactions.  Very  often  the  reaction  rate  Jar 
the  main  reaction  between  the  fuel  and  air  is  so  fast  that  the 
concentration  of  the  main  reacting  species  can  be  calculated 
using  this  assumption  of  local  equilibrium.  Thus,  if  the  reaction 
where  (X)  is  any  third  molecule. 


(02)  + (X)  t 2(0)  + (X) 
are  assumed  in  equilibrium,  then 


(*») 


(0)/{02)1/2  = 4.1  exp  (-29500 /T  )( moles  / cn?  ) ^ ^ (1.) 

Tor  the  temperatures  and  pressures  of  interest  in  many  appli- 
cation.*: the  resulting  value  of  { 0 ) is  very  small  compared  with 
the  concentration  of  the  major  species,  {N^},  (0^)  al1  ^ the  con- 
centration of  fuel  {F}. 


The  pi’oduction  rate  of  NO  as  given  by  (3)  and  (5)  is 


(NO) 


5.74x10 


14 


{N2){(?2} 


1/2 


exp(-67600/T ) 


(G) 


wNQ  = 1 .09.1014p3/2Yt)  Yj/2  exp(-G7600/T)  • (7) 

3 . 3 

in  grams/cm  sec.  Here  p is  the  density  in  gram/cm  and  Y and 


J 


. 
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Yq  are  the  mass  fractions  of  N?  and  0^. 

The  local  chemical  equilibrium  assumption  for  the  reacting 
species  can  be  used  to  calculate  the  local  mass  fraction  o£  0 
in  terms  for  example  of  the  mass  fraction  of  the  elements , "whi eh 
arc  conserved  by  the  chemical  reactions. 

In  previous  analyses  of  the  production  of  110  in  laminar 
and  turbulent  diffusion  flames  the  production  rate  given  above, 
Eq.(7),  was  used  in  connection  with  the  conservation  equations 
for  {NO}  which  were  integrated  numerically  to  provide  the  dis- 
tribution of  {NO}.  Here  we  shall  use  the  fact  that  reaction 
rate  is  very  sensitively  dependent  on  temperature,  because 

the  activation  temperature  T = 67C00°K  of  the  reaction  rate  of 
Eq.(7)  is  very  large  compared  with  T,  to  simplify  the  analysis 
of  "NO"  production  in  laminar  and  turbulent  diffusion  flames. 

We  shall  use  the  thin  flame  assumption  of  Burke  Schumann, 
see  Williams  (1965),  to  calculate  the  distribution  of  temperature 
arid  concentration  of  the  species  N ^ and  0 

II  - "NO"  PRODUCTION  IN  LAMINAR  DIFFUSION  Fl.AMES 

In  figure  1 we  represent  schematically  the  configuration 
of  a typical  diffusion  flame.  A thin  flame,  where  the  main 
reaction  between  the  fuel  and  the  air  takes  place,  separates  a 
region  without  fuel  from  a region  free  from  oxygen.  Also  shown 
in  Fig  1 are  typical  concentration  and  temperature  profiles 
cori’esponding  to  a cross  section  BB'. 

According  to  the  overall  expression  for  the  NO  product  ion 
given  above,  Eq.(7),  this  production  will  only  take  place  at  the 
region  inmediately  adjacent  to  the  thin  flame,  on  the  air  side, 
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whore  the . temperature  has  not  decreased  much  below  the  adiabatic 
flame  value  T^  corresponding  to  the  infinitely  thin  flame.  When 
the  temperature  decreases  slightly  below  T^,  the  NO  production 
becomes  negligible  because  of  its  large  sensitivity  with  tempera- 
ture. We  have  skectched  in  r.ig.l  the  region  of  NO  product -J-on 
adjacent  to  the  thin  flame. 

Because  of  the  small  thickness  of  the  region  of  NO  produc- 
tion we  may  talk  about  the  mass  x^ate  of  NO  produced  per  unit  flame 
surface  mjj  . Then,  we  con  use  this  rate  as  a surface  source  of 
NO,  located  on  the  thin  diffusion  flame,  when  calculating  the  NO 
distribution  in  the  flow  field. 

The  mass  of  NO  produced  per  unit  flame  surface  and  time 

'"no  is  £ivor‘  hy 

mN0  “ ! wNOdy  • ( 8 ) 


where  y .is  the  distance  normal  to  the  thin  flame,  and  the  integral 
is  extended  to  the  region  close  to  the  thin  flame  where  is  sig- 

nificant. To  calculate  m ^ we  need  to  know,  according  to  Eq.(7), 

the  temperature  and  concentration  profiles  f<r  Y,  and  Y 

N 2 ^ 2 * 

(the  pressure  may  be  assumed  to  be  constant  in  low  speed  dif- 
fusion flames).  However,  the  NO  production  will  only  be  signif- 
icant in  a thin  region  close  to  the  thin  flame,  and  therefore 
we  may,  when  evaluating  the  integral  in  Eq.(6),  neglect  the  small 
changes  in  p,  and  Y„  ‘in  the  NO  production  zone,  using  in  Eq.(3) 
their  values  calculated  at  the  thin  flame  (subscript  i).  V.'e  shall 
only  retain  the  variation,  across  the  thin  production  zone,  of 

Yq  and  of  the  Arrhenius  exponential. 

2 

In  the  thin  production  zone  the  mass  fraction  of  0 ^ may 
be  approximated  by 
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. 
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ymo2/pfDf 


(9) 


where  m^  is  the  inass  rate  of  consumption  of  0?  per  unit  flame  sur- 
face and2l)  is  the  diffusion  coefficient  of  0^  in  the  mixture.  The 
Arrhenius  exponential  may  be  linearized  following  I'rank-Kai.ianetsh ii 
around  the  value  T_  because  T /'['  - 67600/T,.  >>  1.  That  is  we 

1 Q i t 

shall  write 

exp(-T/T ) = exp  ( -T  /!'  ) exp{T  (T-T,)/T^}  (10) 

a d j a x i 

where  the  temperature  difference  (T-Tf)  in  the  exponent  will  be 
approximated  in  the  production  region  by 


Tf  = y 


oT 

9y 


-yq 


/kf 


Here  q^  is  the  heat  flux  per  unit  flame  surface  from  the  flame 
toward  ?he  oxygen  side  of  the  flame,  and  k^  is  the  heat  conducti- 
vity of  the  mixture  at  the  thin  flame  conditions. 


As  a result  of  these  approximation  we  obtain,  in  gr/c.n  .sec, 

m"  = 1 . 0 9 p ^ ^ 2 Y f f Y^/^e'^a/idy  ~ 

N0  i,UJpf  Jo  o2 


■So  ■ o.se.io^'hr 


m" 0 1/2  q"0 

(r-“)  ( 


_.i  JjL.) 

kf  T2 


-3/2e-Ta/Tf 


Ol) 


This  result  is  obtained  as  the  first  approbation  when  evaluating 
the  integral  of  Eq.(8)  by  means  of  Laplace  asymptotic  technique 
for  large  values  of  the  exponent  (67600/T^-). 


6 . - 


Taking  into  account  the  fact  that  the  heat  flux  will  in 

most  cases  be  proportional  to  the  oxygen  mast'  consumption  rate  m"  , 

^ 2 

the  NO  production  rate  v/ill  then  be  inversely  proportional  to  m" 
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11 


J-0-  - 0.06.,u’  V’VTf<->  Xt) 


0 


l°2 


0 1/2 
c i 3/2  pfD,-  2 py‘ 

f~. — i)  

m1'  I)  , 

e a 0 f 


T12) 


where  L = p ,-c  D,/k£  is  the  Lewis  number, 
e f p f f 

In  particular,  in  the  analysis  of  the  quasi-steady  burning 

of  fuel  droplets  of  radius  a in  air,  both  mV  and  q"  v.'ill  be 

- 1 U 2 1 2 

proportional  to  a ' . On  the  other  hand,  the  thermodynamic 

properties  at  the  flame  surface  are  found  to  be  independent  of 

the  droplet  radius  when  the  Burke- Schumann , thin  flame,  analysis 

is  used  for  the  main  reaction.  This  analysis  can  be  used  to  predict 

the  flame  characteristics  as  long  an  the  droplet  radius  docs  not 

become  so  small  that  flame  extinction  occurs.  Then 


NO 


-1/2  3/2 
'v  a a = a 


and  the  ratio  between  the  rates  of  NO  production  and 

0^  or  fuel  consumption  is  proportional  to  a,  indicating  the 
importance  of  a reduction  in  droplet  size  in  reducing  the  NO 
product  ion . 


Notice  that  mj^  has  been  found  as  a function  of  the  oxygen 
mass  consumption,  heat  Yelease  rate,  and  other  thermodynamic  data 
of  the  thin  diffusion  flame  corresponding  to  the  main  reaction 
between  the  fuel  and  air.  All  these  data  can  be  obtained  by  using 
the  standard  procedures  of  calculating  diffusion  flames  with  the 
Burke- Schumann  assumption. 
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Ill  - The  production  of  "NO"  in  turbulent  diffusion  flames 

When  analyzing  a turbulent  reacting  flow  we  encounter  the 
problem  that  the  fluctuations  of  the  reacting  species  are  of  the 
order  of  their  mean  values,  so  that  the  mean  production  rate  can 
not  be  written  in  terms  of  the  "Laminar  reaction  rates"  or  rates 
calculated  using  the  mean  values  of  the  concentration  and  tempe- 
rature. This  difficulty  led  Hawthorne,  Weddell  and  Hottel  (1949) 
to  the  introduction  of  a probability  distribution  function  for  th 
concentration  and  a method  of  calculating  the  mean  product  con- 
centration. The  probability  distribution  function  has  been  used 
later  by  Toor  (.1962  ),  O'Brien  (1971),  Fletcher  and  Haywood  (1971) 
Kent  and  Bilgor  (1972  ),  Rhodes  et  al.  ( 1973),  Bus!)  and  Tendelj 
( 1973  ),  I.in  and  O'Brien  (1974)  and  Bilger  and  Kent  ( 1974  ).  See 
also  the  review  by  Bilger  (1976). 

V/e  leave  for  an  Appendix  the  analysis  of  the  turbulent  mixing 
and  infinitely  fast  irreversible  reaction  of  a low  speed  stream 
of  fuel  with  air.  We  show  that,  under  the  assumption  of  equal 
dif f usivities  of  heat  and  mass  , the  instantaneous  concen- 

trations and  temperature  of  all  the  species  can  be  calculated 
in  terms  of  the  instantaneous  concentration  of  one  of  the  species 
i.e.  molecular  oxygen,  where  it  is  different  from  zero. 

Thus  we  obtain  the  relations 


T-T 


20 


1 10+qY10/cp-T20 


T-T20  + q(Y2-Y20>^Cp  ; 


= 0 


(13) 


1 lO-"1  20_qY20/vCp 


Y - Y 
20  1 2 


^ Y M 2 ^ 2 0 


— = 0 


(14  ) 


VY10H  Y 20 


between  the  temperature  T,  the  Nitrogen  mass  fraction  Y and  the 

' 2 


8 


mass  fraction  of  oxygen  Y „ (for  the  notation  see  the  Appendix). 

We  could  similarly  calculate  the  mass  fraction  of  all  the 
other  species  and  thus  the  density  if  the  pressure  fluctuations 
are  neglected  in  the  equation  of  state.  We  could  thus  calculate, 
using  Eq  . ( 7 ) , the  instantaneous  production  rate  of  NO  as  a function 


If  the  probability  distribution  function,  P^tY^),  for  the 
mass  fraction  of  oxygen  were  known  then  we  could  use  the  instan- 
taneous production  rate  w to  calculate  the  local  mean  produc- 


tion rate  a: 

N 0 


'NO 


r 20 

= I u 
1 o 


NO^Y2^P2^Y2^dY2 


(15) 


To  calculate  this  integral  we  shall  make  use  of  the  fact 
that  because  T >>  T the  NO  production  is  so  sensitively  dependent 
on  temperature  that  the  local  NO  production  will  be  limited  to 
those  times  whan  the  concentration  of  Y^  is  sufficiently  small,  so 
that  T is  close  to  the  value  T ^ (adiabatic  flame  temperature) 
obtained  from  Eq.  (13)  by  writing  Y^  = 0.  That  is 


Tf-T 


20 


T10_T20hq^ 10^Cp 


Tf~  r20"qY20/vCp 

1 10_T20_qY20/vCp 


(16) 


2 

when  Y^  is  not  small  so  that  (T^-TjT^/T^  becomes  large  compared 

with  one,  the  Arrhenius  'exponential  in  w,,„  becomes  so  small  that 

If  0 

the  NO  production  may  be  neglected. 

Then  all  the  factors  entering  the  integrand  of  Eq(15)  with 
1 / o 

the  exception  of  |°2I  and  the  Arrhenius  exponent,  can  be 


9 . - 
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taken  as  constant  and  evaluated  at  the  condition  corresponding 
to  Y^  = 0 (and  Y^  = 0)  v/hich  we  shall  designate  with  the  subscript  f. 

Thus  we  obtain 


“no  = 1.09.10i"(p3/2YH^)fF2(0  + )'[ 


(17) 


where 


I = / 20Y^/2  exp ( -T  /T ) dY 


(IP  ) 


To  calculate  the  integral  I we  can  linearize  the  Arrhenius  exponent 
around  the  value  T = Tf,  and  take  into  account  the  relation  (13) 


between  T and  Y ^ written  here  a; 


Tf  - T = (T10“T20  + riYlo/Cp)Y2/(Y20  + vYlo) 


(19) 


so  that  the  .integral  simplifies  to 


-T 


i/Tf  J 


oY2/2  uxp(_ 


^ J~  1 0 1 2 0 * q ^ 10/Cp  ^ Y2 

Y 2 0 + V Y 1 0 


)dY, 


(20  ) 


where  the  upper  limit  in  the  integral  has  been  replaced  by  « , 
because  the  integrand  becomes  extremely  small  when  the  exponent 
becomes  moderately  large  compared  with  one. 


Then 
I = 

T 

a 

1 10~T20+q/i0/cp 

-3/2 

/¥ 

r. xp(-T  /T,  ) 
a i 

(21  ) 

Tf 

Y20+vY10 

2 

We 

have  thus  obtained 

the  local 

mean  "NO"  product 

ion  rate 

in  terms  of  zone  thermodynami c data  the  adiabatic  diffusion  flame 
temperature  T ^ , and  in  terms  of  the  particular  value  P„(0+,x)  of 
the  probability  distribution  function  of  Y . According  to  Eq.(20) 
of  the  appendix 


P2(0 


x) 


P( 


20 


x ) , 


Y20+vY10 


Y20+vY10 


(22) 


so  that  the  NO  production  is  thus  given  in  terms  of  this  particular 
value,  at  Y=f  , of  the  probability  distribution  function  for  an 
inert  species,  of  unit  concentration  in  the  fuel  stream  an  zero 
concentration  in  the  air  stream. 


. The  fact  that  the  production  of  Nitric.  Oxide  is  dependent, 
not  on  smooth  moments,  but  on  the  detailed  chape  of  the  probabi- 
lity distribution  function  for  an  inert  species,  illustrates 
the  need  for  measurements  and  theories  providing  this  p.d.f.. 

The  measurements  of  the  p.d.f.  in  a two-dimensional  turbulent 
mixing  layer  by  Rebc 1 lo ( 1 9 7 2 ) show  that  these  functions  are  not 
Gaussian  but  some -t iines  approach  bimodal  distributions. 
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APPENDIX 


"THE  ROLE  Or  THE  PROBABILITY  DISTRIBUTION  FUNCTION  FOR  THE 
CONCENTRATION  OF  AN  INERT  SPECIES  IN  THE  ANALYSIS  OF  TUR 
BULENT  DIFFUSION  FLAKES". 


We  shall  analyze  in  this  Appendix  some  aspects  of  the 
turbulent  mixing  and  chemical  reaction  between  two  reactants 
undergoing  a one-step  irreversible  reaction  of  the  form 

A1  + n2A2  - n 3 A 3 

We  shall  assume  that  the  molecular  diffusion  velocities 
are  related  to  the  instantaneous  concentration  field  by  Pick's 
law  with  equal  diffusion  coeficients  D for  all  the  species.  We 
shall  assume  that  the  specific  heats  of  all  the  species  are  equal 
to  c , and  constant.  In  addition,  we  shall  limit  the  analysis  to 
l°w  speed  flows,  for  which  the  viscous  dissipation  and  the  work 
associated  with  pressure  gradients  can  be  neglected.  We  shall 
consider  the  thermal  diffusivity  to  be  equal  to  D. 

With  these  assumption;,  the  instantaneous  temperature  and 
concentration  field  will  satisfy  the  following  conservation 
equations 

L (Y.)  = w.  (2) 

t 1 X 

L (c  T 1 
P 

where  L is  the  operator 
L (Y)  = — (pY)  t V.(pvY)  - V.(pDVY) 


w1q 


iE 

at 


(3) 


(4) 


i 


In  these  equations  is  the  mass  fraction  of  species  i, 

. ->• 
p the  density,  T the  temperature,  p the  pressure,  v the  velocity, 

and  is  the  mass  production  rate  of  species  i per  unit  volume. 

If  we  take  into  account  the  overall  stoichiometric  relation 
(1)  we  can  write 

W = w /V  = -w  /(I  + v)  (5) 

where  v = p^lj/M^  is  the  stoichiometric  mass  ratio  oxygen 
(species  2)  to  fuel  (species  1).  M.  is  the  molecular  mass  of 
species  i.  For  an  inert  species  j,  w.  = 0. 

In  Eq.(3),  q is  the  heat  release  per  unit  mass  of  fuel,  or 

q = h°  + vh°  - (1  + v)hg  ( 6 ) 

in  terms  of  the  specific  heats  of  formation  of  the  reacting 
species  and  products. 

Notice  that  the  conservation  equations  (2)  and  (3)  arc 
valid  for  turbulent  flow  as  long  as  we  use  in  them  the  instan- 
taneous values  for  v,  the  thermodynamic  variables  and  w^. 


I • 


There  are  many  reactions  that  can  be  described  by  an 
overall  reaction  of  the  form  (1)  only  if  the  characteristic 
mechanical  time,  in  this  case  given  by  the  smaller  time  scale  of 
the  turbulent  flow,  is  large  compared  with  some  characteristic 
chemical  time  as  required,  for  example,  to  establish  a stationary 
state  in  the  radical  concentrations. 

From  ttie  conservation  equations  (2)  and  (3),  we  deduce, 
if  Eq.(5)  is  taken  into  account  that 


A.  ?. 


i.(y2 

- vY1)  = 0 

(7) 

L(Y3 

+ (1  + v)Y1)  = 0 

CO 

L(T  + 

qYl/cp)  = 0 

(9) 

so  that  the  combinations 

( Y 2 ' vYj),  ty3  +(1  + V)Y1)  and  (T  + qYj/Cp) 

obey  the  same  equation  as  an  inert  species.  For  Eq.(g)  to  be 
valid  the  term  3p/'dt  in  Eq.(3)  should  be  negligible.  This  is  the 
case  for  low  speed  flows  if  the  mean  pressure  outside  the  mixing 
zone  does  .not  vary  with  time. 


For  species  1 we  must  use  an  equation 


L(Y1)  = Wl 


(10) 


where  the  chemical  production  term  must  be  retained.  If,  however, 
we  consider  cases  for  which  the  characteristic  chemical  time  for 
the  reaction  (1)  is  short  compared  with  the  smaller  time  scale  of 
the  turbulent  flow,  we  can  then  substitute  Eq.(10)  by 


Wj  = 0 ( 11  ) 

or  taking  into  account  the  overall  irreversibility  of  the  reaction 
(1),  by 


Y1  Y2=  0 (12) 

showing  that  the  fuel  and  oxygen  do  not  coexist  anywhere  in  the 
flow  field. 


A . 4 . - 


Eq . ( 3 2 ) , 
a conservation 


together  with  the 
equation 


conservation  Eqs.(7) 


to  ( 9 ) , 


L(Y  j ) 


0 


(13) 


for  the  inert  species,  and  the  momentum  equation 


-r“-  + V.(pvv)  = - Vp  + pi  i-V. 
o l m 


I 


and  equation  of  state 


P = pTR(EYi/M.) 


( 14  ) 


(15) 


when  used  with  appropriate  initial  and  boundary  conditions  would 
determine  the  instantaneous  velocity,  temperature  and  concentration 
fields.  because  the  reaction  rates  do  not  appear  in  these  equa- 
tions the  reaction  is  called  diffusion  controlled. 


These  equations  are  to  be  solved  with  appropriate  initial 
and  boundary  conditions.  If  the  resulting  flow  is  turbulent, 
this  is  a formidable  problem.  Our  purpose  here  is  to  show  how 
we  can  obtain  information  for  the  distribution  of  concentrations  an 
temperature  in  diffusion  controlled  combustion  from  details  of 
the  concentration  distribution  for  an  inert  species.  In  addition 
there  are  cases  where  the  heat  release  associated  with  combustion 
will  not  affect  much  the  turbulent  flow.  In  these  cases, appro- 
priate experiments  in  turbulent  mixing  without  combustion  can  be 
used  to  calculate  the  temperature  and  concentration  distributions 

in  diffusion  controlled  combustion. 

* . 

Let  us  consider,  for  example,  the  turbulent  mixing  of  a 
jet  of  fuel  with  a co-flowing  stream  of  air,  as  shown  in  figure 

1, 


a . - 


A. 


and  let  us  assume  that  we  know  the  concentration  Y of  an  inert 
species,  of  unit  concentration  in  the  incoming  stream  of  fuel, 
and  zero  concentration  in  the  air. 


Thus  Y is  a random  function  f (x  , t),  which  wo  shall 
describe  in  terms  of  the  probability  distribution  function  F(Y  , x) 
If  the  overall  process  is  'stationary  the  time  will  not  appear  in 
this  function.  This  function  is  such  that  P(Y  , x)dY  is  the 
probability  of  having  at  the  point  x a concentration  of  the  inert 
species  between  Y and  Y + dY . 


How,  the  combinations 


-V^VY20 

VY10+Y20 


Y3+(1+v)Y1 


and 


( 1 + v ) Y 


10 


T^VVT20 

Tl0+qY10/cp“T20 


Y2+vcpT/q-ucpT2Q/q-Y20 

-Y20+VCpT10/q-VCpT20/q 


(16) 


a . e . - 


satisfy  the  same  transport  equation  and  boundary  conditions  as 

the  inert  species  Y,  so  that  they  have  the  same  probability  dis - 

— ^ 

tribution  function.  That  is,  P ( f j , x)  dfp  is  the  probability  of 

having  the  combination  (vYi  + Y20”Y2  vY10  hY20  ^ take  a Viil«e  be- 
tween fj  and  f^+df^. 


If  we  now  take  into  account  the  equilibrium  relation  (12). 

I 

we  can  now  obtain  the  probability  distribution  function  for  the 
species  Y^,  Y^  and  Y^  and  for  the  temperature  T. 


Thus  to  calculate  the  probability  distributions  function 
Pj(Y^,  x)  for  Y , we  make  use  of  the  fact  that  at  times  when 
Yj  ^ 0 we  have  Y^  = O.and  then  the  combination  f^  given  in  (16) 
is  a function  of  Y1  : = ( vY 1 r Y? 0 ) /( vY 1 Q + Y 20  ) for  whom  we  know 

the  probability  distribution  function  P(fj  , x).  When  Y^  lies 
between  Y^  cind  Y^  + dY  , f^  lies  between  the  corresponding  value 
of  f1  and  f 1+( vdYj ) /( vY10  + Y2Q  ) 


so  that 


vY  + Y 

P (y  -)  = , $) 
VY10+Y20 


VY10+Y20 


(17) 


if  Y ^ o 

when  Y1  = 0 ,,  f*  = ( Y 2Q-Y 2 ) / ( YJ#t vY 1Q  ) < fc 


where 


f 

c 


,20/(I20,''Y10) 


(ie) 


and  the  probability  of  having,  at  a given  point  x,  a value  of  f^ 
lower  than  fc>  that  is  Y^  = 0,  is 


A.  7 


/ C P(Y  , x)dY 


Then  we  may  write  the  probability  distribution  function  P,  for 
Y as 


PjUi  > *>  = 


vY 1 +Y - n f 

F ( — - , x)  + (/0C  P(V  , x)dY]6(Y1)  (13) 


VY10+Y20  vY10+Y50 


j 


Where  the  first  term  gives  the  probability  of  having  a value  of 
Y1  / 0,  and  the  last  term  involving  the-  Dirac  6 function  is  the 
probability  of  having  Y = 0. 

In  a similar  way 


P ( Y 
2 ' 1 2 


x)  = 


Y20+vY10 


P( 


Y 20  ~ Y 2 
Y20+vY 1G 


)+{/J  P(Y  , x)dY}ifY2)  (20) 


To  obtain  the  probability  distribution  function  for  the  temperature 
Pt(T  , x),  we  take  into  account  that  at  times  where  Y ^ = 0 


0 < 


T-T 


20 


T10“T20+qY10/cp 


< f 

c 


with  T ranging  from  1 2Q  to  Tf  • And  at  times  where  Y?  = 0 


1 


» 


\ 


1 > 


T-T2Q-qY20^Cp  > 
T10-T20-qY20/vCP 


f 

c 


with  T ranging  from  , T Q to  Tf 


. 


A . 8 . - 


and 


Tf  = T20  + (T10  " T 2 0 + qY10/cp)fc  ls  the  adiabatic  flame 


temperature.  So  that  P - P,^  + 
with 


( 21 ) 


T-T 


PT2(T  , x)  - 


T10-T20+qY10/cp 


PC- 


20 


T10_T20+qY10/cp 


, x)  (21  a) 


for  0 < (T-T20)/(T10-T20  + qY10/cp)  < f<j 


and  P^  = 0,  outside  the  interval.  (T2q,  T^) 


and 


PTi(T  , x)  = 


T 2 0 T10  + qY20/vCp 


"P  2 0 ~ q Y 2 0 ^ V Cu  -*■ 

PC .12. 20 p___  } x)  (21fc) 


10'  1 20~qY20/vCp 


for 

1 > (T-T20-qY20/vcp)/(T10-T20-qY20/vcp)  > fQ 
and  P^,  = 0,  outside  the  interval  CT^q  , T^) 

In  a similar  way  we  could  write  the  probability  distribution 
function  for  the  product  species  Yg. 

The  probabi lity, distribution  function  for  the  temperature 
could  also  be  obtained  as  follows.  We  first  notice  that  the 
combination 


T + qVVT20 

T10+qY10/cp"T20 


T-T20*q(VY20)/VCP 

T10-T20“qY20/vCp 


a 


(22) 


A.  9.- 


dif fuses 
equation 
a = 0 s o 


like  an  inert  species,  satisfying  the  differential 
L(a)  = 0 with  the  homogeneous  boundary  conditions 
that  the  solution  is  a = 0.  That  is, 


T^VVT20 

T10+qY10/cp"T20 


T-T2otq(Y2-Y 

T10-T20“qY 20 


20)/VCp 


/ v c , 


0 


(23) 


is  a relation  between  the  instantaneous  values  of  T,  Y and 

1 

In  a similar  way  we  may  obtain  the  following  relation 


Y3+(1+v)Y1 


(l-i-v)Y 


10 


Y20-VVY1 

Y20+VY10 


= 0 


( 24  ) 


between  the  instantaneous  values  of  Y , Y and  Y , which  would 

■L  <-  o 

enable  us  to  calculate  the  probability  distribution  functions  P 
and  P2  given  by  (19)  and  (20). 


Also,  if  Yj  is  an  inert  species,  of  zero  concentration  in 

the  fuel  stream  and  concentration  Y_.  in  the  air  stream,  its 

J ° -> 

probability  distribution  function  P.(Y.  , x)  would  be  given  by 


P3 < V1  • *>  ' 


Y . 

1 ° 


P(  1- 


Y . 
1 


Y . 
30 


x) 


(25) 


a result  which  can  also  be  obtained  from  the  fact  that  the 
instantaneous  values  of  Y ^ , Y2  and  Y ^ satisfy  the  relation 


Y . 
._JL 


Y . 

3 ° 


vYltY20-Y2 

VY10+Y2o’ 


= 0 


(26) 


and  using  the  probability  distribution  functions  given  by  (39) 
and  (20). 


A . J 0 . - 


It  is  interesting  to  notice  that  the  relations  (23)  , 

(24)  and  (26)  give  T,  Y and  Y.  in  terms  of  'i'  , ,at  times  when 

J J ^ 

Y^  = 0,  and  similary  in  terms  of  Y ^ at  times  when  Y - o.  The 

mean  concentrations  Y , Y ^ , Yg  and  T can  now  be  obtained  as  a 
. -*  . . 

function  of  x taking  into  account  the  probability  d ir.tr  ib  at  ion 
functions  which  have  been  obtained  above.  Thus 


= /o1C>YlPl(Yl  ’ *)dYl  = /i-{(Y10  + Y20/v)Y-Y20/v}F(Y)dV  (27) 


Y 0 f 

V*)  = 4 Y2P2(Y2’’X-)dY2  = /oC{Y20-Y(Y20  + vT10)}P(Y)dY  (?8) 

! 


The  mean  values  of  Y^  and  T can  be  obtained  more  directly 
from  (23)  and  (24)  which,  obviously,  are  also  valid  for  the 
mean  values  . 


In  conclusion,  if  we  are  able  to  obtain  either  experimen- 
tally or  theoretically  the  probability  distribution  function 
P(Y  , x)for  an  inert  species  of  normalized  concentration  Y,  such 
that  it  is  unity  in  the  fuel  stream  and  zero  in  the  air  stream, 
then,  using  the  ideas  given  above  we  can  calculate  the  probabi- 
lity distribution  function  for  the  reacting  species  as  well  as 
for  the  temperature  if  the  following  main  assumption  can  be  used: 

a)  Equal  thermal  and  mass  diffusivitie  for  el]  the  species 

b)  Infinite  reaction  rate  with  an  irreversible  reaction. 
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